A scalable DG solver for the electroneutral Nernst-Planck equations

The robust, scalable simulation of flowing electrochemical systems is increasingly important due to the synergy between intermittent renewable energy and technologies such as storage chemical manufacturing. high Péclet regime many applications prevents use off-the-shelf discretization methods. In this work, we present a high-order Discontinuous Galerkin scheme for electroneutral Nernst-Planck equations. chosen charge conservation formulation allows specific treatment different physics: upwinding advection migration, interior penalty diffusion ionic species well electric potential. Similarly, enables treatments in preconditioner: AMG potential blocks ILU-based methods advection-dominated concentration blocks. We evaluate convergence rate through numerical tests. Strong scaling results two preconditioning approaches are shown large 3D flow-plate reactor example.